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ABSTRACT 

The phenomenon of cosmic shear, or distortion of images of distant sources 
{Nj ■ unaccompanied by magnification, is an effective way of probing the content and 

state of the foreground Universe, because light rays do not have to pass through 
matter clumps in order to be sheared. It is shown that the delay in the arrival 
^ ■ times between two simultaneously emitted photons that appear to be arriving 

from a pair of images of a strongly lensed cosmological source contains not only 
information about the Hubble constant, but also the long range gravitational 
effect of galactic scale mass clumps located away from the light paths in question. 
This is therefore also a method of detecting shear. Data on time delays among a 
sample of strongly lensed sources can provide crucial information about whether 
extra dynamics beyond gravity and dark energy are responsible for the global 
flatness of space. If the standard KCDM model is correct, there should be a 
large dispersion in the value of H as inferred from the delay data by (the usual 
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procedure of) ignoring the effect of all other mass clumps except the strong lens 
itself. The fact that there has not been any report of a significant deviation from 
the h = 0.7 mark during any of the H determinations by this technique may 
already be pointing to the absence of the random effect discussed here. 



1. Introduction; time delay anisotropy from primordial matter distribution 

ACDM cosmology models the near Universe in terms of (a) the gravity of embedded 
mass clumps in (b) a smooth 'cosmic substratum' of expanding space. While inflation may 
ensure an Euclidean mean geometry, fluctuations caused by the (a) phenomenon operating 
over Hubble scales ought to be observable. It would be very important, therefore, to test if 
the statistical effect of the gravity of virialized structures distributed throughout the near 
Universe exists. Such an effect can manifest itself as a random delay in the arrival times 
of two light signals emitted simultaneously at separate positions and detected by the same 
observer O, after having propagated through different paths. 



A calculation of the (finite) variance in such a delay was provided by Lieu & Mittaz 
(2007), who assumed a smooth Universe perturbed by primordial matter fluctuations of 
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power spectrum P(k). Here we simply sketch the essential steps on how it is done. Denote 
the gravitational perturbation of an otherwise zero curvature Friedmann-Robertson- Walker 
Universe (as inferred from WMAP1 and WMAP3, viz. Bennett et all 2003 and Spergel et 
al 2007) by y), where the x-axis is aligned with the light path and y is a vector along 
some direction transverse to x. If the angle one ray makes w.r.t. the other at O is 9 and 
the comoving light pathlength is D, the relative delay in the arrival conformal time may be 
written as 

r (0) - r (0) = \ I 2x'9-V$(x',0)dx l + \ I x' 2 {9 ■ V) 2 $(x', 0)dx' + ■ ■ ■ , (1) 
c Jo c Jo 

where V is the gradient operator transverse to the vector x, viz. along the y direction. 

The variance in the difference between the delays in the two signals has its lowest order 
term ensuing from the first integral in Eq. (1), as 

\5t{9)} 2 = ^ f° 2x'dx' f D 2x"dx"9 i 9AVMx', 0)V"$(x", 0)), (2) 
c Jo Jo 

where (V^V^) is the correlation function between the two spatial gradients of $, with 
the indices i,j denoting the two transverse directions (0, 1,0) and (0,0, 1), and summation 
over repeated indices is implied. The important point about Eq. (2) is its dependence on 
6> 2 , leading to a standard deviation 5r(9) ~ 9. This zeroth order contribution to 5r, though 
large, has no observational consequence because it depicts a coherent delay 8t{29) = 2St(9) 
which, as will be explained in the next section (see also Bar-kana 1996 and Seljak 1994) 
simply causes a global absolute shift in the angular positions of images without changing 
relative positions. The next order contribution to 5r(9) would come from the second integral 
in Eq. (1), i.e. Sr(9) ~ 9 2 . It depicts the genuinely random excursion in the relative delay 
between the two rays which in principle is observable. 

2. The impossibility of observing the coherent time delay between the light 

curves of strong lensing multiple images 

The best way of demonstrating this impossibility is by means of a concrete example. 
Consider the two-dimensional problem of Figure la, where all light rays are confined to the 
a^-plane. Let a spherically symmetric lensing system be at comoving position (d, 0, 0), with 
the observer at the origin, and let a source S at distance D = 2d cause two images to appear 
on opposite sides of S, say at angular positions 9 + and — 9- (S is slightly off the x-axis). To 
begin with, suppose let were no gravitational perturbations anywhere near the lines of sight. 
Then the distances of closest approach are b± w d9±, assuming that GM <^ b±. Moreover, 
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by the symmetry of the problem, the rays beyond the lens meet the line from lens to source 
at the same angles ± 9±. The angles of deflection of the rays are 4GM/b±. We shall not need 
the actual formulae for position of the source, and the time delay between the two signals. 

Next, in Figure lb we introduce another mass affecting the 'observer's half of the light 
paths, in the plane of both paths, but to one side, bending the rays in the same direction. 
Assume that the mass is not too close, so that its gravitational field may be described by a 
linearly varying potential, of the form $ = —<f)(x)y, where <f)(x) is a smooth function peaked 
in some region of x. Then the null geodesic equation reduces to 

^ = = 1m, 

dx 2 c 2 dy c 2 ' 

with solutions 

1 r 

y± = ± x6± + — / 2(x- x')Mx')dx' . 
c Jo 

The slopes of the two curves when they reach the lens are given by 



y ' ± = ± e± + 2 I (f)(x')dx', 
Jo 



i.e. the angle between them remains at the previous value of 9 + + 0_. We may picture the 
wave fronts as moving backwards from the observer whilst always maintaining orthogonality 
with the direction of propagation. The rays are bent upwards, and the wavefronts in the 
upper ray have less speed than those in the lower one (since the latter experiences weaker 
potential), by just the amount needed to ensure this condition is satisfied. 

Given that what we see in Figure lb is the same as that in Figure la, the lensing mass 
L must also be in a slightly different position, moved upwards along the lensing plane by the 
corresponding amount, 



Sy 



L 



\ f 2(d - x')(f)(x')dx' 
c Jo 



. Unless there are further masses affecting the propagation on the far side of the lens (i.e. 
the parts of the light paths between the lens and the source) the remainder of the diagram 
is exactly as before, except for being rotated by the small angle, 

i> = \ I (P(x')dx' 
c Jo 

as shown in Figure lb. Specifically if the perturbing mass is displaced in the +y (or +z) 
direction, the rotation will be about an axis parallel to z (or y), and in the sense of +y (or 
+z). 
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On the far side, we could similarly trace wavefronts of the signal propagating from the 
source. We can think of the time delay difference as occurring close to L, between wavefronts 
of rays propagating from both ends. Owing to the slight misalignment a between the source 
and lens, this time delay difference will not be zero, but it will be exactly the same for Figure 
lb and Figure la. There is no extra contribution to the difference from the perturbing mass, 
in the context of our lowest order (linear, or coherent) theory. 
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Fig. 1. — The propagation of light signals from a source to an observer, with a strong 
gravitational lens in between. The top sketch depicts the situation where the lens is the only 
disturbance, while in the bottom sketch the effect of a remotely located second mass is also 
present on the plane of the two rays. The dotted lines represent virtual rays going directly 
from source and lens to the observer, with the bottom sketch showing how these rays are 
deflected by the linear perturbation. 
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3. Time delay anisotropy from a statistical ensemble of galaxies 

The analysis in section 1 of light propagation time through a Universe of mass fluctuation 
could be continued with a calculation of the correlation function between potential gradients, 
e.g. Eq. (2) for the lowest order form of such a function, by expressing the integral in terms 
of the matter power spectrum P(k), 

<vj$(0v;$(O> = ilpr / d ^hk^P{k) (3) 

where r = r' — r" and 

with <5$fc being the standard deviation of the potential over length scales 2n/k. 

Now the objective of this paper is to investigate whether time delay measurements can 
probe the mass distribution in the near (z < 1) Universe where non-linear matter clumping 
is important, because similar considerations in the context of primordial matter that fills 
the high redshift Universe have already been made, with the conclusion that such forms of 
matter causes negligible additional delays (Surpi et al 1996, Bar-kana 1996, Seljak 1994). 
One could persist with the approach of Eq. (1) by adopting a modified primordial matter 
spectrum P(k) that includes an 'extension' to the non-linear, or large k, regime, (Peacock 
& Dodds 1994, Smith et al 2003), except that concerning the effect of widely spaced and 
compact clumps on two closely separated light rays it actually makes more sense to calculate 
directly the time delay induced by the gravitational perturbation of a random ensemble of 
clumps. This is because of several reasons: (a) while the correlation function of Eq. (3) 
is relatively simple for the unobservable lowest order term of 5t (i.e. sections 1 and 2) 
the observable next order term has a much more complex form; (b) the shape of P(k) is 
unreliable on sub-Mpc scales; (c) the use of a Poisson clump distribution is justified by the 
'nearest neighbour' interaction phenomenon. More elaborately the differential time delay 
between two closely spaced rays is due mostly to proximity galaxies located at distances 
sufficiently small from the two rays in question for mass clustering (or compensation) to 
exert any significant modification - see below. 

Now even if mass correlation over larger scales can be neglected by appealing to Poisson 
clumps, there still is a form of P(k) to represent this kind of matter inhomogeneity (Eq. 23 
of Peebles 1974). Owing to reason (a) above, however, it is much easier to start from first 
principles. Let us return to the one-way Shapiro time delay when light skirts a mass clump 
m at 'impact parameter' y, or more precisely transverse position y (with components along 
the conventional (0,1,0) and (0,0,1) directions, both perpendicular to i = x) w.r.t. the 



7 



unperturbed ray, which is 



2Gm 



In 




(5) 



where x s and xi s are the observer-source and clump-source distance respectively. Take a pair 
of rays with separation 8y <C y, the difference in delay between them is 



If geometry is globally flat, the comoving separation Sy between the two rays at any position 
of comoving distance x from the observer O is given by 



where 6 is the angle subtended at O by two point sources of distance d apart, both being 
at the same comoving length x s away from O, to which our two rays map back. These two 
points could even mark the positions of a pair of strong lensing images, when the light rays 
associated with the images are perturbed by external mass clumps, in which D,, 
the distance to the lensing plane, and 5r is the relative delay in the light arrival times 
between the two images. 

The lowest order observable effect is the incoherent contribution to the relative delay 
between the two rays from each clump, which originates from the second spatial derivative 
of the clump's gravitational potential. From section 2 and Eq. (1) we that this is the 
contribution arising from the 5r ~ 2 ~ {Sy) 2 terms of Eq. (6). In this light, it is clear that 
only the 2nd and 3rd terms on the right side of Eq. (6) are relevant. Thus, when we square 
the equation to form the variance, we obtain 



where the angle averages employed to go from Eq. (6) to Eq. (8) were (cos 2 1?) = 1/2 and 
(cos 4, #) = 3/8, with d being the angle between y and Sy. 

For the rest of this section we shall indeed focus our attention upon one manifestation 
of shear: the delay in photon arrival times between two strong lensing images. Under this 
scenario a random walk arises as a result of the two rays skirting all the other clumps on 
either side of the light path. Their loci then become like two long snakes (Hamana et al 
2005, Gunn 1967a,b) as the rays are deflected, largely in tandem, though there is always a 



St = 





(7) 
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small and random relative change in directions, which is the same shear phenomenon as the 
incoherent relative delay caused by each clump. Both these relative processes of deflection 
and delay accumulate like from one clump to the next if the clump distribution is 
Poisson, as assumed. 

Thus, specifically concerning time delay, as the rays continue their journey the variances 
(^ r )hicoherent fr° m & U the clump encounters add, so that to evaluate the total excursion in 
the arrival time difference one must perform a cylindrical integration with the axis along the 
x-direction. There is however one subtlety here. Since y is the comoving impact parameter 
its value for deflections at finite z was smaller by by the factor 1 + z. Fortunately this factor 
cancels out when all the components of the integrand are assembled, though the same does 
not happen when one computes image distortions by shear, as we shall see. The cumulative 
variance for a Poisson ensemble of clumps of comoving number density n (i.e. neglecting 
with caveats the evolution of clump properties, see below) is obtained by an integration down 
the light path to be 

,„,., (Gm\ 2 f D ' c , f y ™2irydy 2tt / GmV nDf . 
[5r mcohcrcnt (9)] 2 = 2 (^-3-) J o {Syfndx jf -JUL = - -^-9\ (9) 

if one considers only the contribution to {St) 2 from foreground mass clumps that interact 
with the light rays during redshifts z < z\. In arriving at the final expression of Eq. (9) use 
was made of Eq. (7). We may get rid of the dependence on y min by assuming that it equals 
the comoving distance at which one or more clumps satisfy the condition y < y m i n , i.e. 

2/min = 7=f- (10) 

The contribution from background (zi < z < z s ) clumps may likewise be calculated and 
included. One would then arrive at a total variance of 

c 2 (5r incohcrcnt ) 2 = ^ (^p) Q 2 cl Df(Df + Dl)6\ (11) 

after employing the relation Gnm = J2 i Gn i m i = 3HqQ c \/ (8ii), with Q d being the mass 
density of clumps as a fraction of the critical density. 

We intend to pursue an application of the above development, by predicting the effect 
of external field galaxies on the time delay between strong lensing images, for comparison 
with observations. Before doing so, however, several remarks about (Sr) 2 andom are in order. 
Apart from the most obvious fact that its final form scales only with one property of the 
clumps, viz. Q c \, Eq. (11) is valid in the limit 5y y. In terms of the comoving separation 
d between sources, Eq. (7), this implies (since 5y < d and y min < y) that d <^ y m i n . Now the 
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images to be used as testbeds have a typical comoving separation d < 50 kpc, i.e. a scenario 
under which the criterion is satisfied, because 50 kpc is not much larger than the size of a 
galaxy. More elaborately, if y min < d one would 'see' a galaxy when looking at the sky along 
any direction: we assume this is not the case. 

The next remark is that the sole role played by proximity clumps can be seen from 
the dependence of (6r) 2 on y min and not y max . Thus, on the question of relative time delay 
caused by galaxies one does not need to take account of large (Mpc scale or more) distances 
over which mass compensation by galaxy clustering is important. This justifies a posteriori 
our use of a random clump ensemble without appealing to P(k). It also provides the reason 
why the effect of voids on the time delay fluctuations can be ignored: since the maximum 
separation between the two rays is small compared with the typical inter-clump spacing 
the void-to-void accumulation of the randomly varying second spatial derivative of the void 
potential function is completely negligible over such distance scales (over much larger (CMB 
acoustic) distance scales this phenomenon could play a role in lensing deflections via mass 
clustering, see Holz & Wald 1998, Seljak 1996). The inclusion of void effects will in principle 
add further signal to the variance 5r inco herent for the presently assumed (and justified) Poisson 
mass distribution, Thus even though the neglected contribution is small, it does mean that 
our estimate of 5r incohercnt is conservative. 

4. Testing the dynamics of global geometry by strong lensing time delay 

The interpretation of cosmological time delay data is usually confined to considerations 
of the delay within the strong lens system and its immediate environs, with the overall aim 
of inferring the Hubble constant H from the observations (Refsdal 1964). If e.g. the lensing 
mass distribution is a singular isothermal sphere, the relative time delay between two images 
at angular distances 9 a and 9 b from, and on opposite sides of, the axis of symmetry is given 



In Eq. (12) it is assumed, of course, that A and B are images of the same source, usually 
a time variable background quasar. The Hubble constant clearly affects the delay via the 
distance dependence, viz. Atab ~ DiDs/Dis ~ Hq 1 (other cosmological parameters also 
play a role because D is a multi-dimensional function, but their effects are minor, as noted by 
Grogin & Narayan 1996). Hence time delay measurements via light curve alignment between 
images A and B of the quasar, coupled with knowledge of redshifts, can in principle lead to 
a determination of H . 



by 




(12) 
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To date approximately ten strong lensing systems with time delay measurements are 
available (Saha et al 2006), in each case the delay between two images separated by several 
arcseconds is typically found to lie within the 10 - 100 days range. One pair of such multiply 
lensed quasars with similar parameters, SDSS J1004+4112 and HE0435-1223, were reported 
by Fohlmeister et al 2006 and Kochanek et al 2006 respectively. We shall employ this pair to 
illustrate how the cosmological distribution of galaxies near the light path can substantially 
enlarge the random uncertainty in the value of H . For SDSS J1004+4112 where the smallest 
observed image separation, between images A and B, was 9 ~ 4 arcsec and the comoving 
distances are D s = 4.77 Gpc, Di = 2.45 Gpc in an Vl m = 0.3, = 0.7 cosmology, Eq. 
(11) gives <5Ti nco h crG nt ~ 24 days if we persist with the ACDM standard model by adopting 
its breakdown of the matter budget to set Q c \ = 0.15 (specifically this assumes that half the 
baryons, hence approximately the same fraction for dark matter also, of the low z Universe 
resides in galaxies and their halos, see Fukugita 2004 and Fukugita et al 1998). In fact, 
taking the above parameters as typical, one could proceed to recast Eq. (11) into a more 
convenient form: 



^Tineoherent = 29. 1 ( h/0.7) 2 (ft d /0. 15) (0/5 arcsec) 2 ( A/2.5 Gpc) f 



1 + 



A 



days. (13) 



Since the observed image delay of 38.4 ± 2.0 days is on par with 5r ran dom, any estimation 
of Ho that attributes all the observed delay to physics within the strong lens system would 
have caused this value to vary randomly from the truth by almost 100 %. Other images of 
SDSS J1004+4112 can also be used as testbeds: the separation of image B from C (also A 
from D) is 9 pa 20 arcsec, while the expected time delay between them is > 560 days (> 800 
days for A and D), due solely to the strong lens system itself. From Eq. (13) we see that 
once again 5ri n coherent is comparable to these delays because of its 9 2 scaling. 

A repetition of the above analysis to HE0435-1223, where D s = 6.44 Gpc, A = 1.74 
Gpc and 9 pa 2 arcsec, results in a similar though less drastic conclusion, viz. 5ri nco herent ~ 
4.5 days versus the observed delay of 14.4 ± 0.8 days. The random error in Ho here should 
then account for an additional fluctuation 8H /H pa 30 %. 

In summary, the distribution, evolution, and mass budget of galaxies as understood 
in the context of the standard cosmological model leads to the prediction of a random (or 
incoherent) relative delay between the light arrival times from two images of a strongly lensed 
background quasar comparable with the observed delay. Since the latter has routinely been 
interpreted as an effect caused principally by the gravitational field of the lens, and moreover 
a value of H is derivable from it if perturbations outside the strong lens are absent, the 
question of whether additional and hitherto unknown dynamics are responsible for the global 
flatness of space could be addressed by examining the statistical variation in the H values 
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that emerge from a large number of strong lensing delay measurements, when such a database 
becomes available. 

If this variation distributes around a value of H that agrees with other methods of 
determination, with a standard deviation matching the expectation from Eq. (10), it would 
imply that for the first time the ensemble gravitational effect of many galaxies spread over 
cosmological distances has been detected. If, on the other hand, the variation distributes 
tightly around the accepted value of H with no room for extra perturbations, the possibility 
of a new physical phenomenon that complements (even replaces) the law of gravity as the 
distance scale becomes large must then be inevitable. The fact that there has not been any 
report of a significant deviation from the h = 0.7 mark in the value of H as determined by 
this type of time delay technique may already be pointing to the absence of the random effect 
discussed here. If this turns out to be really the case, the stability of the global geometry as 
revealed here would present a serious challenge to cosmology. 

We end by pointing out that the possible influence of foreground matter on the mea- 
surement of H was considered in a recent work (Fassnacht et al 2006) under the scenario 
of this matter being clumped into several foreground groups of galaxies. The present paper, 
however, calculates for the first time the effect of a random ensemble of many foreground 
clumps on strong lensing time delay; we then demonstrated that this introduces a large scat- 
ter in the resulting value of H . The cause of such a scatter stems mainly from light skirting 
clumps without passing through them, i.e. cosmological time delay data contain precious 
information on shear. 

The author thanks an anonymous referee for very helpful criticisms towards improving 
this paper. 
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